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This paper presents a procedure for decomposing an nth order filter 
into cascadable second order sections. The procedure is optimal in that it 
minimizes the maximal response range for the sections within the frequency 
band of interest. The procedure, based on a modified version of the Bottle- 
neck Assignment Algorithm, describes methods of listing all the optimal 
decompositions as well as of finding a special "nested" optimal decom- 
position. 

I. INTRODUCTION 

Let <f)(s) be a transfer function 

♦« = % 

where / and g are polynomials with real coefficients, and the degree of 

/ ^ degree of g. 

We consider all the decompositions of the form tf>(s)=0i(s)0 a (s) • • • 

<f} t (s) where 

*» - m m 

/,(«) and <7,(s) are real polynomials and the degree of /,• does not exceed 
the degree of g t . The gi are quadratic polynomials, except when the 
degree of g is odd; then one g t is linear. 

Let L be a passband region for 4>, where L is a finite union of passband 
intervals. Then for every <t> { , a number d(<f>i) is defined by 



<2(*<) = 20 log.o 



Max | 4>i{jw) 
Min I tj>i(ju) 
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Also let 

d = Max d(<f>i). (3) 



Then d is a function of the decomposition. We present a procedure 
that determines the decomposition (s) with a minimal d. E. Lueder 
proposed this optimality criterion. 1 

II. METHOD 

First, we artificially equate the number of zeros [zeros of /(s)], and 
poles [zeros of g(s)], by adding a suitable number of "zeros at infinity" 
corresponding to constant unit polynomials. Next we make this mutual 
number even by adding a zero and a pole at infinity, if necessary. In 
this way we get, say, 2t zeros and 2t poles. 

Pairing two zeros creates an /, ; a real zero can be paired with any 
other real zero, while a complex zero must be paired with its conjugate 
in order to get a real /,• . The same is true for creation of g t by pairing 
of poles. 

In the following we assume that all poles, except perhaps one, are 
complex and therefore fixed paired. We call the real zeros which are 
not fixed paired free zeros. 

Next we make all possible pairings of the free zeros. Each such 
pairing, together with the fixed pairing, decomposes /(s) and g(s): 

m -«•)/.« ••■/.(«); 

g(s) = g 1 (s)g 2 (s) •■■ g t (s). 
Then we compute the matrix D = (d ik ), where the elements 



, - M (4) 



d 

are computed from definition (2). The element d ik represents the "cost" 
of matching zero-pair i with pole-pair k. 

An assignment is a feasible set of matchings. Using the Bottleneck 
Assignment Algorithm, we determine an assignment ki , • • • , k t for 
which 

Max d iki 

i - 1 t 

will be minimal. We call this minimum the optimal d value for this 
pairing of free zeros. Going through all the possible pairings of free 
zeros, we find an optimal pairing which yields the smallest optimal 
d value. 

Since an optimal assignment (for a given optimal pairing) is usually 
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not unique, procedures for obtaining all the optimal solutions (assign- 
ments) or a nested solution are given. A nested solution is obtained by 
taking an optimal solution, fixing the matching with the largest d 
value, and then proceeding to look for an optimal assignment for the 
remaining t— 1 //s and t— 1 g k 's, and so on. 

III. THE BOTTLENECK ASSIGNMENT ALGORITHM 

This section discusses the Bottleneck Assignment Algorithm and its 
adaptation to the present problem. 

Let U = (w,,) be a real t X t matrix. A matching is an ordered pair 
of integers (i k , j k ) l^ik^t, l^j k ^t. We associate with the matching 
iik , jk) the corresponding cell in U. The element in this cell u it is 
called the cost of the matching. 

A set A = {(»* , j k ); k = 1, • ■ • , t] of t matchings (cells) is called an 
assignment if in every row and in every column of U there is a cell that 
belongs to A. The bottleneck assignment problem looks for an assign- 
ment which minimizes the maximum of its matchings' costs. 

The Gross algorithm 2 , is based on the following iterative step: 

(*) An assignment A and a real number a, which does not exceed 
all the costs of A, are given; then either a new assignment A' is con- 




A=ARBITRARY ASSIGNMENT 



n=MAXIMAL COST OF 
MATCHINGS OF A 



APPLY THE ITERATIVE STEP (*) 




NO 



A=A' 



Fig. 1 — Flow chart for solving the bottleneck assignment problem. 
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structed, so that a exceeds more costs in A' than it does in A, or it is 
established that no such A' exists. 

The flow chart in Fig. 1 solves the bottleneck assignment problem. 
This algorithm is fast and requires small memory space, since only the 
present assignment must be stored. 
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Fig. 2 — Flow chart for finding the optimal pairing and its optimal d-value. 
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The basic algorithm was modified to find the optimal pairing of the 
free zeros and its optimal d-value (Fig. 2). Note that the value of a in 
every iterative step (*) does not exceed the minimum of the optimal d 
values of the checked pairings. Thus any pairing that does not reduce 
the d value already obtained is immediately disregarded. 

Also, for each pairing we use the optimal assignment of the preceding 
pairing, as an initial assignment. Thus the costs of the initial assignment 
matchings that correspond to the fixed paired zeros do not exceed the 
current d value. These procedures considerably reduce the amount of 
computation required for finding the optimal pairing and its optimal d 
value. 

IV. CREATION OF THE NESTED SOLUTION 

Let U denote the cost matrix which corresponds to the optimal 
pairing. The nested solution is created by successively applying the 
bottleneck assignment algorithm t times and modifying U each time 
in such a way that the matching with the largest cost becomes fixed 
and irrelevant in the further computations. 

Let (i% , j* k ) be the matching with the largest cost at a certain stage. 
Then (i* k , j* k ) becomes a part of the nested solution. U is then modified 
as follows: 

U it \, = co for all s ^ ft ', 

Ut,i t * = °° for all s 7* i* ; 

U it .. it . = 0. 

It is easy to verify that this modification has the properties described. 

V. A COMPUTATIONAL METHOD TO GENERATE ALL THE OPTIMAL 

ASSIGNMENTS 

Let U denote again the cost matrix which corresponds to the optimal 
pairing, and let d* be the optimal d value. We call a cell (i, j) admissible 
if u tt ^ d*. The problem of listing all the optimal assignments then 
becomes the problem of listing all possible assignments that use only 
admissible cells. Using the flow chart of Fig. 3 can accomplish this. 
The number of operations can be seen to be dependent on the order of 
the columns of U. The dependence is quite complicated. However, a 
good rule of thumb for reducing the number of operations is to re- 
arrange the columns in ascending order according to the number of 
their admissible cells. 
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I(l) = 0,I(2)=0 f ...,I(t) = 



I(j) = o 



j=2 



j = j - 1 




PRINT OUT THE ASSIGNMENT 
[(10), 0.(1(21,2), (iCtt.t)] 



Fig. 3 — Flow chart for generating all optimal solutions. 
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